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ABSTRACT 

This paper investigates the flow and heat mass transfer past an infinite vertical plate with variable thermal 
conductivity, heat source and chemical reaction. The non-linear, coupled partial differential equations together 
with the boundary conditions are reduced to dimensionless form. The resulting equations are discritized using 
implicit finite difference scheme of Crank-Nicolson type and solved numerically. The velocity, temperature and 
concentration profiles are presented graphically with tabular presentations of the Skin friction, rate of heat and 
mass transfer which are all computed and discussed for different values of parameters of the problem. 
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I. INTRODUCTION 

The study of effects of porous boundaries on flow and heat transfer with mass transfer is important 
because of many engineering design applications in the field of chemical and geophysical sciences. Permeable 
porous plates are used in the filtration processes and also for a heated body to keep its temperature constant and 
to make the heat in solution of the surface more effective. The study of stellar structure on solar surface is 
connected with the mass transfer phenomena. Its origin is attributed to difference in temperature caused by the 
non-homogeneous production of heat, which in many cases can rest not only in the formation of convective 
currents but also in violent explosions, [1]. 

Most of the practical situations demand for fluids that are non-Newtonian in nature which are mainly 
used in many industrial and engineering applications. It is well known that a number of fluids such as molten 
plastic, polymeric liquid, food stuffs, etc, exhibit non-Newtonian character. The boundary layer flow of non- 
Newtonian fluids over a starching sheet has been studied extensively in the recent years. [2] presented a work on 
flow and heat transfer in power law fluid over a stretching porous surface with variable surface temperature. 
Very recently, [3] considered the effects of buoyancy and variable thermal conductivity in a power law fluid past 
a vertical stretching sheet in the presence of non-uniform heat source. 

The study of convective fluid flow with mass transfer along a vertical porous plate in the presence of 
magnetic field and internal heat generation is receiving considerable attention due to its useful applications in 
different branches of Science and Technology such as cosmical and geophysical sciences, fire engineering, 
combustion modeling, etc. [4] analyzed the MHD free convective flow past an accelerated vertical porous plate 
by finite difference method. 

[5] investigated the transient free convection in cold water past an infinite vertical porous plate. [6] 
discussed free convection and mass transfer flow through a porous medium past an infinite vertical porous plate 
with time dependent temperature and concentration while [7] investigated the effect of chemical and thermal 
diffusion with Hall current on unsteady hydromagnetic flow near an infinite vertical porous plate and [8] 
reported the effects of applied magnetic field on transient convective flow in a vertical channel. [9] analyzed the 
effect of Hall current MHD free convection flow along an accelerated porous heated plate with mass transfer and 
internal heat generation. [10] studied transient free convection flow of a viscous dissipative fluid past a semi- 
infinite vertical plate. [11] investigated the unsteady MHD convective heat transfer past a semi-infinite vertical 
porous moving plate with variable suction. Also, [12] analyzed the effect of mass transfer in unsteady MHD flow 
and heat transfer past an infinite porous vertical moving plate while [13] discussed the unsteady free convective 
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flow and mass transfer of a rotating elastico - viscous liquid through porous media past a vertical porous plate 
and [14] studied the unsteady free convection and mass transfer boundary layer flow past an accelerated infinite 
vertical porous plate with suction. [15] investigated unsteady MHD free convective flow and heat transfer along 
a vertical porous plate with variable suction and internal heat generation. [16] obtained the analytical and 
numerical results for free convection flow along a porous plate with variable suction in porous medium while 
[17] obtained the numerical approach on parameters of the thermal radiation interaction with convection in 
boundary layer flow at a vertical plate with variable suction and [18] presented the convective heat and mass 
transfer flow over a vertical plate with Nth order chemical reaction in a porous medium. 

The aim of the paper is to investigate the unsteady heat mass transfer flow past an infinite vertical plate 
with variable thermal conductivity and heat source in the presence of chemical reaction. The governing 
equations are solved numerically using the implicit finite difference scheme of Crank-Nicolson type. The effect 
of the parameters on the velocity, temperature and the concentration distributions of the flow fields are discussed 
and shown through graphs while the skin friction, rate of heat and mass transfer are discussed and shown using 
tables. 



II. PROBLEM FORMULATION 

We consider an unsteady infinite vertical isothermal porous plate of laminar natural convection flow of 
dissipative and radiating fluid in the presence of transverse magnetic field surrounded on one side by infinite 

mass of fluid like air or water and both at same temperature T' Sj and the mass concentration C x initially. At time 
t'> 0, the plate temperature and the mass concentration is raised to T and C , causing the presence of 

temperature and concentration difference T — T x and C w — C x respectively. As the plate is infinite in extent, 

the physical variables are functions of y' and t' where y' is taken normal to the plate and the x'-axis is taken along 
the plate in the vertically upward direction. It is assumed that both the variable thermal conductivity and the nth 
order chemical reaction are not constant. Under the Boussinesq's approximation, the governing equations in this 
case [19] and [20] for the flow are continuity, momentum, mass concentration and energy respectively (which is 
an extension of [21]: 

dv' n 

— = 0 (1) 
By' 

dt' dy' dy' 2 V xJ 



8T' ,dT Id 
+ v 



dt' dy' pC dy' 
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The initial and boundary conditions relevant to the fluid flow are: 
t'<0, u' = 0,T' = Tl,C' = Cl for ally' 

t'>0, u' = 0,T' = T; v ,C' = C' w aty' = 0 (5) 

u -^0,T' -^Z,C -»C1 asy'^oo 
where V is the kinematic viscosity of the grey fluid, <7 is the Stefan-Boltzmann constant, B 0 is the constant 
magnetic field intensity, p is density, K is the permeability, g is the gravitational constant, [3 is the thermal 
expansion coefficient, /?* is the concentration expansion coefficient, T' is the temperature, C" is the mass 
concentration, D is the chemical molecular diffusivity, R is the chemical reaction, C p is the specific heat at 
constant pressure, q r is the radiative heat flux, u and v' are velocity components in x and y directions 
respectively, t is the time, K(T' ) is the variable thermal conductivity, b y is the Forchheimer parameter of the 
medium, b is the Hall current parameter, Q is the volumetric rate of heat generation and n is the reaction order 
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while T is the wall temperature, is the free stream temperature, C' w is the species concentration at the plate 

surface and C' x is the free stream concentration. 

From the continuity equation (1), it is clear that the suction velocity is either a constant or a function of time. 
Hence, on integrating equation (1), the suction velocity normal to the plate is assumed in the form, 

V = — V 0 where V 0 is a scale of suction velocity which is non-zero positive constant. The negative sign 

indicates that the suction is towards the plate and V 0 > 0 corresponds to steady suction velocity normal at the 
surface. 

Assuming the radiative heat flux from the Rosseland approximation to have the form: 

^ = -w(C-r 4 ) (6 ) 

dy v ' 

o is the Stefan-Boltzmann constant, a * is the mean absorption effect for thermal radiation constant. We assume 
that the temperature differences within the flow are sufficiently small such that T' 4 can be expanded in a Taylor 
series about and neglecting higher order terms give: 

T' 4 = 4T'Z 3 -3T' 4 (7) 

The variable thermal conductivity depends on temperature. It is used by [22], as follows: 

K(T') = k 0 {l + r(T'-T:)} (8) 

where k 0 is the thermal conductivity of the ambient fluid and y is a constant. 

III. METHOD OF SOLUTION 

To solve the governing equations in dimensionless form, we introduce the following non-dimensional quantities: 



u= jL,, = yOx, t= 



u 



0 



U 



V 



C'-C 



(T'-T'Y 



u- 



Sc = — ,Pr = -^,Ec = 

d K c p (t;-t:) 



Gr= — m — ' Gc= — tz? — ' a= w: {9) 



K = 



K'U n 



-,N = 



u 



k Q U 0 



,M = 



PU 0 2 



u- 



- b s 



Qq 2 

k Q U Q 



The governing equations on using (9) into (1), (2), (3), (5) and (6) to (9) into (4) reduce to the following: 



T^^govergigg eqg$JjSns usingl(9^ into (1), (2), (3), (5) and using (6) to (9) into (4) reduce to the following 
a = — r- Mh \U + Gr6 + GcC-bJJ ( 
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subject to the boundary conditions 

t<0, U =0,6 = 0,C = 0 
t>0, U =0,6 = l,C = l 
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(10) 



(11) 



(12) 



(13) 



as y — >oo 

where Pr is the Prandtl number, Sc is Schmidt number, Ec is Eckert number, Gr is thermal Grashof number, Gc 
is mass Grashof number, bj is the inertia number, and M is magnetic field, K is porosity, N is the radiation, a is 
suction, T is the variable thermal conductivity, Kr is the chemical reaction, b is Hall current and S is heat source 
parameters. 



www.theijes.com 



The UES 



Page 79 



Heat Mass Transfer Flow past an Infinite Vertical Plate. 



IV. NUMERICAL PROCEDURE 

In order to solve the unsteady, non-linear, coupled equations (10) to (12) under the 
boundary condition (13), an implicit finite difference scheme of the Crank-Nicolson type has 
been employed. The finite difference equations corresponding to equations (10) to (13) are as 
follows: 

-rjj£ + (1 + 2i; pr l - rp£ = r 2 UU + (l-2r a - yr 3 - r 4 )Uj + (r 2 + yr 3 )U ] M + AtGr0 i i 

2 ( 14 ) 

+MGcCI - Mb x [U j) 
-r^ +(Sc + 2r 1 )cr 1 - r x C£ = tfiU + (Sc - 2r 2 - yr 3 Sc)C/ +(r 2 + yr 3 Sc)q +i 

/ \n (15) 

-KrScteyq) 

-qr$* + (Pr+ 2qr x ) 6>/ +1 - qrfi£ = qr 2 0/_ x + (Pr- 2qr 2 - yr 3 Pr- NAt + SAt)0/ 

2 2 2 (16) 
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The mesh sizes along y- direction and time t-direction are Ay and At respectively 
while the index i refers to space y and j refers to time t . The finite difference equations 
(14) - (16) at every internal nodal point on a particular n-level constitute a tridiagonal system 
of equations which are solved by using the Thomas Algorithm. 

In each time step, the concentration and temperature profiles have been computed first from 
equations (15) and (16) and then the computed values are used to obtain the velocity profile 
which meets the convergence criteria. 

The skin friction coefficient, rate of heat and mass transfer in terms of Nusselt number and 
Sherwood number respectively are given by 

C f = 



(su} 


, Nu = 




, Sh = 
















y=0 


\ d yj 


y=0 




y=0 



V. RESULTS AND DISCUSSION 

The numerical solutions are simulated for different values of the Prandtl number (Pr), Schmidt number 
(Sc), Eckert number (Ec), thermal Grashof number (Gr), mass Grashof number (Gc), Inertia number (bi), 
radiation (N), magnetic field (M), porosity (K), variable thermal conductivity ( X ), suction ( a ), reaction order 
(n), chemical reaction (Kr), Hall current (b) and heat source (S) parameters. The following parameters values are 
fixed throughout the calculations except where otherwise stated, Pr = 0.71, Sc = 0.62, Ec = 0.01, M = 1.0, K = 
1.0, Gr= 1.0, Gc = 1.0, N = 0.1,1)! =1.0, b = 1.0, S = 1.0, a =1.0, X = 0.1, Kr = 0.1, n= 1.0. 

The velocity profiles are illustrated in Figures 1 to 12 for different values of Prandtl number (Pr = 0.71, 
1.0, 3.0, 7.0), Schmidt number (Sc = 0.22, 0.62, 0.78, 2.63), thermal Grashof number (Gr = 1, 5, 10, 15), mass 
Grashof number (Gc = 1, 5, 10, 15), Inertia number (b, = 1, 10, 20, 30), magnetic field (M = 1, 5, 10, 15), 
porosity (K = 0.1, 0.5, 1.0, 1.5), radiation (N = 1, 5, 10, 15), suction ( a = 2, 4, 6, 8), variable thermal 
conductivity (T = 0.1, 0.5, 1.0, 1.5), chemical reaction (Kr = 0.1, 1.0, 10, 100) and heat source (S = 3, 5, 7, 9) 
parameters. 
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In Figure 1, it is seen that the velocity decrease with increasing Prandtl number also Figure 2 reveals 
that, the velocity decrease with increase in the Schmidt number and likewise, Figure 3 shows that, the velocity 
decreases with increase in the magnetic field parameter. But, Figure 4 indicates that, the velocity increase 
whenever porosity parameter increases also Figure 5 represents the thermal Grashof number and it is observed 
that, the velocity increase with increasing thermal Grashof number just like Figure 6 which reveals that, the 
velocity increase when the mass Grashof number increased. While from Figure 7, it shows that, the velocity 
decrease with increasing radiation parameter and similarly for the suction parameter, Figure 8 indicates that, the 
velocity decrease with increase in the suction parameter. In Figure 9, it is seen that, the velocity increase with 
increasing variable thermal conductivity parameter whereas in Figure 10 we found that, the velocity decrease 
when the chemical reaction parameter increase just as in Figure 1 1 where the velocity decreases with increasing 
Inertia number. Figure 12 indicates that, the velocity increase whenever the heat source is increased. 

The concentration profiles are illustrated in Figures 13 to 15 for different values of Schmidt number (Sc 
= 0.22, 0.62, 0.78, 2.63), suction (a = 2, 4, 6, 8) and chemical reaction (Kr = 10, 30, 50, 70) parameter as in 
Figures 13, 14, and 15 respectively. 

In Figure 13, it is noticed that, the concentration decrease with increase in the Schmidt number and 
likewise Figure 14 reveals that, the concentration decrease with increase in the suction parameter. Similarly, 
Figure 15 indicates that, the concentration decrease whenever the chemical reaction parameter increased. 

The temperature profiles have been studied and presented in Figures 16 to 20 for different values of 
Prandtl number (Pr = 0.71, 1.0, 3.0, 7.0), radiation (N = 1, 5, 10, 15), suction ( a = 2, 4, 6, 8), variable thermal 
conductivity (T = 0.1, 0.5, 1.0, 1.5) and heat source (S = 3, 5, 7, 9) parameters shown in Figures 16, 17, 18, 19 
and 20 respectively. 

In Figure 16, it is observed that, the temperature decrease with increase in the Prandtl number. Also, 
Figure 17 reveals that the temperature decrease with increase of the radiation parameter and similarly for the 
suction parameter, Figure 18 indicates that, the temperature decrease whenever the suction parameter increased. 
Figure 19 shows that, the temperature increase with increase of the variable thermal conductivity parameter and 
similarly for Figure 20, the temperature increase with respect to increase in the heat source parameter. 

Tables 1 to 3 are the tables for Skin friction, Nusselt number and Sherwood number for the numerical 

solution. 

The Skin friction is illustrated in Table 1 for different values of Prandtl number Pr, thermal Grashof 
number Gr, mass Grashof number Gc, Schmidt number, chemical reaction Kr, suction a , magnetic field M, 
radiation N and S heat source parameters. In the Table, it is observed that increase in the Prandtl number, 
Schmidt number, chemical reaction, suction, magnetic field and radiation parameters leads to decrease in the 
Skin friction while the Skin friction increases whenever the thermal Grashof number, mass Grashof number and 
heat source parameter increases. 

The Nusselt number is presented in Table 2 for different values of Prandtl number Pr, thermal Grashof 
number Gr, radiation N, Hall current b and a suction parameters. From the Table, it is noticed that the Nusselt 
number increase with increasing Prandtl number, radiation and suction parameters but increasing the thermal 
Grashof number and Hall current parameter results to decrease of the Nusselt number. 

Similarly, the Sherwood number is demonstrated in Table 3 for different values of Schmidt number Sc, 
chemical reaction Kr and a suction parameters. Results from the Table shows that the Sherwood number 
increases whenever the Schmidt number, chemical reaction and suction parameters increase. 
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Figure 1 : Variation of Velocity against y for different values of Pr. 
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Figure 2: Variation of Velocity against y for different values of Sc. 




y 



Figure 3: Variation of Velocity against y for different values of M. 
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Figure 4: Variation of Velocity against y for different values of K. 
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Figure 5: Variation of Velocity against y for different values of Gr. 




y 

Figure 6: Variation of Velocity against y for different values of Gc. 
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Figure 7: Variation of Velocity against y for different values of N. 
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Figure 11: Variation of Velocity against y for different values of bi. 
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Figure 12: Variation of Velocity against y for different values of S. 
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Figure 13: Variation of Concentration against y for different values of Sc. 
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Figure 14: Variation of Concentration against y for different values of a . 
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Figure 15: Variation of Concentration against y for different values of Kr. 
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Figure 16: Variation of Temperature against y for different values of Pr. 
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Figure 17: Variation of Temperature against y for different values of N. 




Figure 18: Variation of Temperature against y for different values of a . 
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Figure 19: Variation of Temperature against y for different values of X . 
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Figure 20: Variation of Temperature against y for different values of S. 



TABLE 1: Skin friction for different values of Pr, Gr, Gc, Sc, Kr, a , M, N and S. 



Pr 


Gr 


Gc 


Sc 


Kr 


a 


M 


N 


s 


Q 


0.71 






0.62 


0.1 






0.1 




0.6439 


1 






0.62 


0.1 






0.1 




0.6058 


0.71 






0.62 


0.1 






0.1 




1.9552 


0.71 






0.62 


0.1 






0.1 




1.8857 


0.71 






0.78 


0.1 






0.1 




0.6227 


0.71 






0.62 


1 






0.1 




0.6320 


0.71 






0.62 


0.1 






0.1 




0.5407 


0.71 






0.62 


0.1 






0.1 




0.5399 


0.71 






0.62 


0.1 






1 




0.6249 


0.71 






0.62 


0.1 






0.1 


3 


0.6979 



TABLE 2: Nusselt number for different values of Pr, Gr, N, b and CC . 



Pr 


Gr 


N 


b 


a 


Nu 


0.71 


1 


0.1 


1 




0.8356 


1 


1 


0.1 


1 




1.1964 


0.71 


5 


0.1 


1 




0.8180 


0.71 


1 


1 


1 




1.1653 


0.71 


1 


0.1 


10 




0.8163 


0.71 


1 


0.1 


1 


4 


2.2608 



TABLE 3: Sherwood number for different values of Sc, Kr and CC . 



Sc 


Kr 


a 


Sh 


0.62 


0.1 


l 


1.1604 


0.78 


0.1 


l 


1.3499 


0.62 


1 


l 


1.3640 


0.62 


0.1 


4 


2.4522 



VI. CONCLUSION 

The paper presents the heat mass transfer flow past an infinite vertical pate with variable thermal 
conductivity, heat source and chemical reaction. The dimensionless governing equations are non- 
dimensionalized and then solved numerically using the implicit finite difference scheme of Crank-Nicolson type. 
Numerical solutions are presented for the fluid flow and heat mass transfer characteristics for different values of 
parameters involved in the problem. The present study will serve as a scientific tool for understanding more 
complex flow problems concerning with the various physical parameters. 
The conclusion shows that: 

❖ The velocity increase with increase in the thermal Grashof number, mass Grashof number, variable thermal 
conductivity, porosity and heat source parameters whereas the velocity decrease with increasing Prandtl number, 
Schmidt number, Inertia number, magnetic field, radiation, suction and chemical reaction parameters. 
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♦> The temperature increase with increase in variable thermal conductivity and heat source parameters while the 

temperature decrease with increasing Prandtl number, radiation and suction parameters. 

♦♦♦ The concentration decrease with increase in Schmidt number, suction and chemical reaction parameters. 
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